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Abstract 

In this paper, we consider a generalized second order nonlinear ordinary differential equation of 
the form x + (k\x q + £2)2; + k^x 2q+l + kiX q+1 + Xix = 0, where fcj's, i = 1,2,3,4, Ai and q are 
arbitrary parameters, which includes several physically important nonlinear oscillators such as the 
simple harmonic oscillator, anharmonic oscillator, force-free Helmholtz oscillator, force-free Duffing 
and Duffing-van der Pol oscillators, modified Emden type equation and its hierarchy, generalized 
Duffing-van der Pol oscillator equation hierarchy and so on and investigate the integrability prop- 
erties of this rather general equation. We identify several new integrable cases for arbitrary value 
of the exponent q, q £ R. The q = 1 and q = 2 cases are analyzed in detail and the results are 
generalized to arbitrary q. Our results show that many classical integrable nonlinear oscillators can 
be derived as sub-cases of our results and significantly enlarge the list of integrable equations that 
exist in the contemporary literature. To explore the above underlying results we use the recently 
introduced generalized extended Prelle-Singer procedure applicable to second order ODEs. As an 
added advantage of the method we not only identify integrable regimes but also construct integrat- 
ing factors, integrals of motion and general solutions for the integrable cases, wherever possible, 
and bring out the mathematical structures associated with each of the integrable cases. 
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I. INTRODUCTION 



A. Overview of the problem 

In a recent paper 1 we have shown that the force-free Duffing- van der Pol (DVP) oscillator, 

x + (a + (3x 2 )x - 7X + x 3 = 0, (1) 

is integrable for the parametric restriction a = -5 and 7 = — -J5. In Eq. (JI} over dot 
denotes differentiation with respect to t and a, (3 and 7 are arbitrary parameters. Under 
the transformation 

w = —xe~p\ z = e~p\ (2) 
Eq. (Q) with restriction a = 4 and 7 = — J? was shown to be transformable to the form 

w" - ^w 2 w' = 0, (3) 

which can then be integrated 1 . 

In a parallel direction, while performing the invariance analysis of a similar kind of prob- 
lem, we find that not only the Eq. ((TJ but also its generalized version, 
4 3 

x + (— + [3x 2 )x + —x + x 3 + 5x 5 = 0, <5 = arbitrary parameter, (4) 

p p 

is invariant under the same set of Lie point symmetries 2 . As a consequence one can use the 
same transformation (J2J to integrate the Eq. Q. The transformation (J2J) modifies Eq. (@J 
to the form 

«/' - y«) 2 «)' + 5W 5 = (5) 

which is not so simple to integrate straightforwardly. However, we observe that this equa- 
tion coincides with the second equation in the so called modified Emden equation (MEE) 
hierarchy, investigated by Feix et al. 3 , 

x + x l x + gx 2l+1 = 0, / = l,2,...,n, (6) 

where g is an arbitrary parameter. 

In fact Feix et al. 3 have shown that through a direct transformation to a third order 
equation the above Eq. (JHJ) can be integrated to obtain the general solution for the specific 



choice of the parameter g, namely, for g = „ + L 2 . For this choice of g, the general solution 
of (0) can be written as 

, , / (2 + 3l + l 2 )(t + h) 1 \* T T 
x(t) = { Kt + Ii y +1 + {2 + 3l + m ) . h, h : arbitrary constants. (7) 

Consequently Eq. (jlj) can be integrated under the specific parametric choice 8 — ^, and 
it belongs to the / = 2 case of the MEE hierarchy (jHJ) with g = Now the question arises 
as to whether there exist other new integrable second order nonlinear differential equations 
which are linear in x and containing fifth and other powers of nonlinearity. As far as our 
knowledge goes only few equations in this class have been shown to be integrable. For 
example, Smith 4 had investigated a class of nonlinear equations coming under the category 

x + f(x)x + g{x) = 0, (8) 

with f(x) = (n + 2)bx n — 2a and g(x) = x(c + (bx n — a) 2 ) where a, b, c and n are arbitrary 
parameters. He had shown that the Eq. (jHJ) with this specific forms of / and g admits explicit 
oscillatory solutions. However, one can also expect that there should be a number of inte- 
grable equations which also admit solutions which are both oscillatory and non-oscillatory 
types in the class 

x + (kix q + k 2 )x + k 3 x 2q+1 + k 4 x q+1 + X x x = 0, q G R, (9) 

where fcj's, i — 1,2,3,4 and Ai are arbitrary parameters. When q — 1, Eq. becomes the 
generalized MEE 

'x + {kix + k 2 )x + k 3 x 3 + k^x 2 + X x x = 0, (10) 
and for q = 2 it becomes 

x + (hx 2 + k 2 )x + k 3 x 5 + k 4 x 3 + X t x = 0. (11) 

We note that Eq. (jlj) is a special case of (JTTJ). 

Needless to say Eq. © is a unified model for several ground breaking physical systems 
which includes simple harmonic oscillator, anharmonic oscillator, force-free Helmholtz oscil- 
lator, force-free Duffing oscillator, MEE hierarchy, generalized DVP hierarchy and so on. 

As noted earlier there exists no rigorous mathematical analysis in the literature for the 
second order nonlinear differential equations which contain fifth or higher degree nonlinearity 

3 



in x and linear in x and the results are very scarce on integrability or exact solutions. Our 
motivation to analyze this problem is not only to explore new integrable cases/systems of 
Eq. but also to synthesize all earlier results under one approach. 

Having described the problem and motivation now we can start analyzing the integra- 
bility properties of Eq. JSJ). To identify the integrable regimes we employ the recently in- 
troduced extended Prelle-Singer procedure applicable to second order ODEs 5 ~ n . Through 
this method we not only identify integrable regimes but also construct integrating factors, 
integrals of motion and general solution for the integrable cases, wherever possible. 



B. Results 



We unearth several new integrable equations for any real value of the exponent q in 
Eq. Q. In the following we summarize the results for the case q = arbitrary only and 
discuss in detail the q = 1, q = 2 and q = arbitrary cases separately in the following 
sections. 

For the choice q = arbitrary we find that the following equations are completely inte- 
grable (after suitable reparametrizations), all of which appear to be new to the literature: 

x + {hx q + {q + 2)k 2 )x + hhx^ 1 + {q + \)k\x = (12) 

x + {(q + 2)k x x q + k 2 )x + klx 2q+l + hk 2 x q+1 + X lX = (13) 

x + (q + 4)k 2 x + k 4 x q+1 + 2{q + 2)k 2 2 x = (14) 
x + {{q + l)k lX q + k 2 )x + ^~^[{q + l)k\x 2q+1 

+ (q + 2)k 1 k 2 x q+1 + k 2 2 x] = 0, r^O (15) 
x + ((q + l)k lX q + (q + 2)k 2 )x + (q + l)[^^lkjx 2q+1 

+hk 2 x q+1 + k\x\ =0, r ^ 0, (16) 

where ki, k 2 , k 4 , Ai and r are arbitrary parameters. We stress that the above results are 
true for any arbitrary values of q. We discuss the special cases, namely, q = 1 and q = 2 
separately in detail in sections 3 and 4 in order to put the results of q arbitrary case in 
proper perspective. 

We show that the Eq. (|12j) is nothing but a generalization of the Dufnng-van der Pol 
oscillator Eq. (0). In a recent work 1 ' 9 three of the present authors have established the 
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integrability of Eq. (|12|) with q = 2. However, in this work we show that the generalized 
Eq. (JT2*)l itself is integrable. Eq. (jl3j) is nothing but the generalized MEE among which the 
hierarchy of Eq. (jHJ), studied by Feix et al. 3 , can be identified as a sub-case. In fact the 
general solution constructed by Feix et al., Eq. (J7J), can be derived straightforwardly as a 
sub-case. Eq. (|13|) also contains the family of equations studied by Smith 4 . In particular 
the latter author have derived general solution for the case k\ < 4Ai, which turns out to be 
an oscillatory one. However, in this work, we show that even for arbitrary values of &2 and 
Ai one can construct the general solution. Interestingly, the system (fTlj) generalizes several 
physically important nonlinear oscillators. For example, in the case q = 1 and 2, Eq. (jl4)l 
provides us the force-free Helmholtz and Duffing oscillators, respectively, whose nonlinear 
dynamics is well documented in the literature 12-16 . Here, we present certain integrable 
generalizations of these nonlinear oscillators. Eq. (j!5|) admits a conservative Hamiltonian 
for all values of the parameters r, k\ and &2 and any integer value of q. We also provide 
the explicit form of the Hamiltonian for all values of q. As a result we conclude that it is a 
Liouville integrable system. As far as Eq. (fTT)|) is concerned we construct a time dependent 
integral of motion and transform the latter to time independent Hamiltonian one and thereby 
ensuring its Liouville integrability. 

The plan of the paper is us follows. In the following section we briefly describe the ex- 
tended Prelle-Singer procedure applicable to second order ODEs. In Sec. IHIl we consider 
the case q = 1 in Q and identify the integrable parametric choices of this equation through 
the extended PS procedure. To do so first we identify the integrable cases where the system 
admits time independent integrals and construct explicit conservative Hamiltonians for the 
respective parametric choices. We then identify the cases which admit explicit time depen- 
dent integrals of motion. To establish the complete integrability of these cases we use our 
own procedure and transform the time dependent integrals of motion into time independent 
integrals of motion and integrate the latter and derive the general solution. In Sec. IIV[ 
we repeat the procedure for the case q = 2 in Eq. (jHJ) and identify the integrable systems. 
In Sec. we consider the case q = arbitrary in (jHJ) and unearth several new integrable 
equations and their associated mathematical structures. Finally, we present our conclusions 
in Sec. EH 
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II. GENERALIZED EXTENDED PRELLE-SINGER (PS) PROCEDURE 



In this section we briefly recall the generalized extended or modified PS procedure before 
applying it to the specific problem in hand. Sometime ago, Prelle and Singer 5 have pro- 
posed a procedure for solving first order ODEs that admit solutions in terms of elementary 
functions if such solutions exist. The attractiveness of the PS method is that if the given 
system of first order ODEs has a solution in terms of elementary functions then the method 
guarantees that this solution will be found. Very recently Duarte et al. 7 ' 8 have modified the 
technique developed by Prelle and Singer 5 ' 6 and applied it to second order ODEs. Their 
approach was based on the conjecture that if an elementary solution exists for the given 
second order ODE then there exists at least one elementary first integral I(t, x, x) whose 
derivatives are all rational functions of t, x and x. For a class of systems these authors have 
deduced first integrals and in some cases for the first time through their procedure 7 . Re- 
cently the present authors have generalized the theory of Duarte et al. 7 in different directions 
and shown that for the second order ODEs one can isolate even two independent integrals 
of motion 9-11 and obtain general solutions explicitly without any integration. This theory 
has also been illustrated for a class of problems 1,9-11 . The authors have also generalized 
the theory successfully to higher order ODEs 10,17 . For example, in the case of third order 
ODEs the theory has been appropriately generalized to yield three independent integrals of 
motion unambiguously so that the general solution follows immediately from these integrals 
of motion 17 . 

We stress that the PS procedure has many advantages over other methods. To name 
a few, we cite: (1) For a given problem if the solution exists it has been conjectured that 
the PS method guarantees to provide first integrals. (2) The PS method not only gives the 
first integrals but also the underlying integrating factors, that is, multiplying the equation 
with these functions we can rewrite the equation as a perfect differentiable function which 
upon integration gives the first integrals directly. (3) The PS method can be used to solve 
nonlinear as well as linear second order ODEs. (4) As the PS method is based on the 
equations of motion rather than on Lagrangian or Hamiltonian description, the analysis is 
applicable to deal with both Hamiltonian and non-Hamiltonian systems. 
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A. PS method 



Let us rewrite Eq. (^J) in the form 

x = -((k lX q + k 2 )x + k 3 x 2q+1 + k 4 x q+1 + X lX ) = <j>(x, x). (17) 

Further, we assume that the ODE (jT7j) admits a first integral I(t, x, x) = C, with C constant 
on the solutions, so that the total differential becomes 

dl = I t dt + I x dx + I x dx = 0, (18) 

where each subscript denotes partial differentiation with respect to that variable. Rewriting 
Eq. (fTTjl in the form (fidt — dx = and adding a null term S(t, x, x)x dt — S(t, x, x)dx to the 
latter, we obtain that on the solutions the 1-form 

(^) + Sx^jdt - Sdx - dx = 0. (19) 

Hence, on the solutions, the 1-forms (|18|) and ()19p must be proportional. Multiplying f)19|) 
by the factor R(t, x, x) which acts as the integrating factors for l)19|) . we have on the solutions 
that 

dl = R((f> + Sx)dt - RSdx - Rdx = 0. (20) 

Comparing Eq. (fTHj) with (j2Uj) we have, on the solutions, the relations 

I t = R(<f> + xS), I X = -RS, h = -R. (21) 

Then the compatibility conditions, It x = I x t, h± = I±t, Ix± = I±x-, between the Eqs. (|21j) . 
provide us 

S t + xS x + cj>S A = -<t> x + <p x S + S 2 , (22) 

Rt + xR x + 4>R X = -{4> x + S)R, (23) 

R x - SR X - RS X = 0. (24) 

Solving Eqs. (}22|) - (j24"j) one can obtain expressions for S and R. It may be noted that any 
set of special solutions (S, R) is sufficient for our purpose. Once these forms are determined 
the integral of motion I(t, x, x) can be deduced from the relation 



fx — r 2 — 



R + Tx {n - r2) 



dx, (25) 



where 

ri = J R(<f> + xS)dt, r 2 = J{RS+—r 1 )dx. 

Equation (|23|l can be derived straightforwardly by integrating the Eq. (J2TJ). 

The crux of the problem lies in finding the explicit solutions satisfying all the three 
determining Eqs. (|22 |) -(|24 |) . since once a particular solution is known the integral of motion 
can be readily constructed. The difficulties in constructing admissible set of solutions (S, R) 
satisfying all the three Eqs. (j22j) - (j21|) and possible ways of obtaining the solutions have been 
discussed in detail in Ref. 9. 

III. APPLICATION OF PS PROCEDURE TO EQ. TO 

Let us first consider the case q = 1 in Eq. @ or equivalently (fTUj) 

x + (kix + k 2 )x + k 3 x 3 + k^x 2 + Xix = 0. (JTUJ) 

Eq. (fTU|) itself includes several physically important models. For example, choosing fcj = 
0, i = 1, ...4, we get the simple harmonic oscillator equation and the choice kx, k 2 = gives 
us the anharmonic oscillator equation. When kx, k^ = Eq. (JTUJ) becomes the force-free 
Duffing oscillator equation 12 . The choice k 2 , ki, X± = provides us the MEE 18 . In the limit 
ks = k 4 = Eq. (jlUJl becomes MEE with linear term which is another linearizable 
equation which we have studied extensively in Refs. 9 and 19. The restriction kx, k 3 = 
leads us to the force-free Helmhotz oscillator 12,13 . In the following we investigate whether 
the system (I10JI admits any other integrable case besides the above. 

We solve Eq. (jlOj) through the extended PS procedure in the following way. For a given 
second order equation, (jlUj) . the first integral I should be either a time independent or time 
dependent one. In the former case, it is a conservative system and we have I t = and in 
the later case we have J ( ^ 0. So let us first consider the case It = and determine the null 
forms and the corresponding integrating factors and from these we construct the integrals 
of motion and then we do extend the analysis for the case I t ^ 0. 



S 



A. The case I t = 



1. Null forms 

In this case one can easily fix the null form S from the first equation in (|21j) as 

g _ -0 _ {{kix + k 2 )x + k 3 x 3 + k A x 2 + Aiz) 
i x 



2. Integrating Factors 

Substituting this form of S, given in into (j2SJ) we get 

+ — ((kix + k 2 )x + k 3 x 3 + k^x 2 + Xix)R± 

= ({k lX + k 2 ) + ^ lX+ k 2 )± + k^ + k A ^+x lX ) ^ R (2?) 

Equation (|2"7)l is a first order linear partial differential equation with variable coefficients. 
As we noted earlier any particular solution is sufficient to construct an integral of motion 
(along with the function S). To seek a particular solution for R one can make a suitable 
ansatz instead of looking for the general solution. We assume R to be of the form, 

R= (A(x) + B(x)xY' (28) 

where A and B are functions of their arguments, and r is a constant which are all to be 
determined. We demand the above form of ansatz, (|28jh due to the following reason. To 
deduce the first integral I we assume a rational form for /, that is, / = ^jjfy ; where / and g 
are arbitrary functions of x and x and are independent of t. Since we already assumed that 
/ is independent of t, we have, I x = ^ x9 ~J 9x and I± = f ±9 ~/ 9± . From (|2*T|) one can see that 
r = j ± = Ua-fSi s = L, = f*9-fgx anc j RS = I X} so that the denominator of the function S 

9 l± Jx9 J9x 

should be the numerator of the function R. Since the denominater of S is x (vide Eq. (J2HJ)) 
we fixed the numerator of R as x. To seek a suitable function in the denominator initially 
one can consider an arbitrary form R = h ^±) • However, it is difficult to proceed with this 
choice of h. So let us assume that h(x, x) is a function which is polynomial in x. To begin 
with let us consider the case where h is linear in x, that is, h = A(x) + B(x)x. Since R 
is in rational form while taking differentiation or integration the form of the denominator 
remains same but the power of the denominator decreases or increases by a unit order from 



that of the initial one. So instead of considering h to be of the form h = A(x) + B(x)x, 
one may consider a more general form h = (A(x) + B(x)x) r , where r is a constant to be 
determined. Such a generalized form of h and so R leads to several new integrable cases as 
we see below. 

Substituting (|2*8|) into ()27|) and solving the resultant equations, we arrive at the relation 

r(x(A x + B x x) + (f>B) = (A + Bx)(f> x . (29) 

Solving Eq. (jSH)) with (ft = — ((kiX + k^x + k^x^ + k^ + Xix), we find nontrivial forms for the 
functions A and B for two choices, namely, (i) k\, k 2 arbitrary and (ii) k\ = arbitrary, k 2 = 
with restrictions on other parameters as given below. The respective forms of the functions 
and the restriction on the parameters are 

(i) ki, k 2 : arbitrary 

(r — l)bo ,k\ o 

A(x) = ^~2 X ~^ 2X )> B\x) = b = constant, r = constant, 

6o(r-l) 36 (r - 1) & (r - l)fc 2 2 
fc 3= 2r2 k u h = — 2 hk 2 , Ai = ^ , (30a) 

(iz) fci = arbitrary, k 2 = 

2r Ki 

k 3 = — k\, k A = 0, Ai = arbitrary parameter (here). (30b) 

2r 2 

We note that the case (ii) cannot be derived from case (i) by taking k 2 = 0. For example, 
choosing k 2 = in (|30aj) we get not only /c 4 = but also Ai = whereas in the case (u) we 
have the freedom Ai = arbitrary, so the cases ()30a|) and (j30b|) are to be treated as separate. 
Making use of the forms of A and B from Eqs. ()30a|) and (j30b|) into (|28|h the integrating 
factor, 'i?', for the two cases can be obtained as 

(i) k\, k 2 : arbitrary 



r V 2 

(w) fci = arbitrary, k 2 = 



r ^ (31a) 



/?= TY , r^O. (31b) 
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We note that bo is a common parameter in the above and it is absorbed in the definition of 
l R\ see Eqs. (J23J) and While deriving the above forms of R (Eqs. (jSTal) and (I31bl) ) we 
assumed that r ^ and for the choice r = we obtain consistent solution only if both the 
parameters k± and k 2 become zero. Of course, this sub-case can be treated as a trivial one 
since when ki,k 2 = the damping term in Eq. (|1(J|) vanishes and the system becomes an 
integrable anharmonic oscillator. In this trivial case we have the integrating factor of the 
form: 

(Hi) ki, k 2 = 

R = x, r = 0. (31c) 



Finally one has to check the compatibility of forms S and R with the third Eq. (|24p . We 
indeed verified that the sets 

B ((hx + k^x + ^^ + ^x' + k^x)) 

W s = 1 > 

X 

R = — , k\, hi = arbitrary, r ^0 (32a) 

{^{^ + k 2X )+±y 

(kixx + ^-r-kjx 3 + Xix) 

(it) S = , 

x 

X 

R = , . — -, ki = arbitrary, k 2 = 0, r^O (32b) 

( + TT + 

and 

/...n o (hx 3 + k 4 x 2 + Xix) . . . 

(m) S = , R = x, ki, k 2 = 0, (32c) 

x 

satisfy the Eq. (|2~%jl individually. As a consequence all the three pairs form compatible sets 
of solution for the Eqs. (22j)-(j24j). 



3. Integrals of motion 



Having determined the explicit forms of S and R one can proceed to construct integrals 
of motion using the expressions (|25)1. The parametric restrictions ()3()aj) and (J30bj) fix the 
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equation of motion (|10|) to the following specific forms, 



x + (kix + k 2 )x + 



(r-1) 



2r 2 



A; 2 x 3 + 3^1^^ + 2k l 2 x 1=0, r ^ 



(in) 



(r-l)fc? 
272 



x 3 + Aix = 0, r^O, 



x + A; 3 x 3 + k^x 2 + Xix = 0, r = 0. 



(33a) 

(33b) 
(33c) 



In the above ki, k 2 , k^, k±, Ai and r are arbitrary parameters. 

We note that the transformation x = y — |p transforms equation ()33a|) to the form 



y + hyy + 



(r-1)*?, (r-l)k 2 



2r 2 



-y - 



2r 2 



y = 0, r^O. (34) 
. So, hereafter, we consider 



2r 2 



Eq. (J31|) is obtained from Eq. (|33bJ) by fixing Ai = 
Eq. ()33a|) as a special case of Eq. (|33bj) and so discuss only Eq. (|33bjl as the general one. 



It may be noted that Eq. (j33b|) includes several known integrable cases. For example, the 
choice r = 3 and Ai = in Eq. (j33b|) yields the MEE 18 . On the other hand the choice 
r = — 1 leads us to the equation x + k\xx — kfx 3 + X\x = which can be solved in terms of 
Weierstrass elliptic function 20 . The other choices of r lead to new integrable cases as we see 
below. 

Substituting the forms of S and R (vide Eqs. (|32b|) and (|32c)l ) into the general form of 
the integral of motion ()25|) and evaluating the resultant integrals, we obtain the following 
time independent first integrals for the cases ()33b|) and (j33cj) : 



[iia] 



h 



x + 



(r-1 
2r 



-kix 2 + 

h 

r 2 X, 



■ 1 ■ 1 2 

x(x + -x + {r _^ h . 



+ (r-l)A 2 + 



r 2 A x 



(r- l)fci 



(n&) 



h 



Akix 



kix 2 + 4fcii + 8A1 



log(A;{x 2 + 4A;ii + 8Ai), r = 2, 



ii = x H — -£ 2 — — ^ log(fcii; + Ai 



2" ^' 

/ 1 = y + ^ 4 + |x 3 + ^ 2 , r = 0. 



(35a) 
r^0,l,2, 

(35b) 
(35c) 
(35d) 



Note that in Eq. (|35ajl . r can take any real value, except 0, 1,2. In the above integrals I\ 
given by Eqs. ()35a|) - 
Eq. l|5Bcj) . 



cj) correspond to the ODE ()33b|) . while ()35d|) corresponds to the 
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Due to the fact that the integrals of motion (|35p are time independent, one can look for 
a Hamiltonian description for the respective equations of motion. In fact, we obtain the 
explicit Hamiltonian forms for all the above cases. 



4- Hamiltonian Description of it!5\) 

Assuming the existence of a Hamiltonian 

I(x, x) = H(x, p) = px — L(x, x), (36) 

where L( X j X ^ IS the Lagrangian and p is the canonically conjugate momentum, we have 

dl dH dp . dL dp . 

dx dx dx ^ dx dx 

m = dH = dp i _dL 

dx dx dx dx 



From (|37j) we identify 



p = I -7-dx, 



X 



Plugging the expressions (jHEJ) into (J38|) one can evaluate the canonically conjugate mo- 
mentum and the associated Lagrangian as well as the Hamiltonian. They read as follows: 



(a) The canonical momenta : 

(iia,b) p = — j— (x+ ^~ ^ —x 2 + ^] , r ^ 0, 1 (39a) 
r — 1 \ r 2 k\ ) 

[tie] p = log(A;ii; + Ai), r = l (39b) 
(Hi) p = x, r = 0. (39c) 

(Note in the above r = 2 is included in Eq. (j39b|) itself). 

(b) The Lagrangian : 

(r — 1) kix 2 rAi x 2 



(*•> L - j2^w^) \ i+ r^r + t) • r ^ ' 1 ' 2 (40a) 

(Hb) L = log(4A;ii: + 8A X + k 2 x 2 ), r = 2 (40b) 

(He) L = ^-\og(k 1 x + Xi) + x(\og(kix + Xi) — 1) — ]-kix 2 , r = l (40c) 

(ill) L = y-y* 4 -y:r 3 -yX 2 , r = 0. (40d) 



13 



(c) The Hamiltonian : 



(iia) H 



r — l)p 



r-2 
r-1 



(r-2) 



(r - 1) 2 rAi, 



2r 



fei 



(mo) ii = — — p H a; p + log( 

k\ 4 p 



(iic) // = -^(e p -A lP +^ 2 - Ax), r = l 

/■•••\ rr P 2 , ^3 4 , ^4 3 , Ai 2 

zm if = x H x H x , r = 0. 

v ; 2 4 3 2 ' 



0,1,2 (41a) 
(41b) 
(41c) 
(41d) 



One can check that the Hamilton's equations of motion are indeed equivalent to the appro- 
priate equation ([TUj) . 

and ()33cj) admit time independent Hamiltonians they can be classified as 



Since Eqs. 

Liouville integrable systems. The important fact we want to stress here is that for arbitrary 
values of r, including fractional values, the equation \33b\l is integrable. 



5. Canonical transformation for the Hamiltonian Eqs. 141]) 

Interestingly, we also identifed suitable canonical transformation to standard particle in 
a potential description for the Hamiltonians (|4Tj) . Now introducing the canonical transfor- 
mations 



x 



x 



2rP 
felt/) 



fcl ■ 



V 



.Mil 

Ar 



, r^0,l, 



p = —k±U, r = 1 



(42) 
(43) 



the Hamiltonian H in Eq. (|41j) can be recast in the standard form (after rescaling) 



lp2 , (1-0 / (r-l)fci^ 
2^ * r (r-2) \ 4r 



(■—2) 
r-l 



+ 



(r-l)Ai TT 2 



\P 2 + ±U 2 + \og yTll . 
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# = < 



_p2 + e - kl U + ^TJ^ 



k 2 



ip 2 + k±TjA + ^ C/ 2 ; 



f/ 2 , r^0,l,2 



r = 2 



r = 1 



r = 0. 



(44) 
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It is straightforward to check that when U and P are canonical so do x and p (and vice 
versa) and the corresponding equations of motion turn out to be 

(2-r) 

fir — Dfci \ ^7 (3-t-) (r — 1) Ai 
U-2( { 4r ) ^ + 2 ' r^0,l (45a) 

U - he' 17 + k x M = 0, r = 1 (45b) 

t> + A; 3 C/ 3 + AiC/ = 0, r = 0. (45c) 

One may note that the equations of motion now become standard type anharmonic oscillator 
equations. 



B. The case I t ± 

In the previous sub-section we considered the case I t = 0. As a consequence S turns out 
to be However in the case It = 0, the function S* has to be determined from Eq. f!22L 
that is, 

St + xS x — ({k\x + k 2 )x + k 3 x 3 + k^x 2 + Xix)S± 

= (hx + 3k 3 x 2 + 2k 4 x + Ai) - (hx + k 2 )S + S 2 . (46) 

Since it is too difficult to solve Eq. (f4T)j) for its general solution, we seek a particular solution 
for 5", which is sufficient for our purpose. In particular, we seek a simple rational expression 
for S in the form 

s = afax) + b{t,x)x 
c(t, x) + d(t, x)x' 

where a, b, c and d are arbitrary functions of t and x which are to be determined. Of 
course, the analysis of this form alone does not exhaust all possible cases of interest. We 
hope to make a more exhaustive study of Eq. (|4^j) separately. Substituting (}4Tj) into ({HI}) 
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and equating the coefficients of different powers of x to zero, we get 
db x — bd x — kid 2 = 0, 

db t — bd t + cb x — bc x + a x d — ad x — 2k 1 cd — (3k 3 x 2 + 2k 4 x + Xi)d 2 

+(kix + k 2 )bd - b 2 = 0, 

cb t — bc t + da t — ad t + ca x — ac x — k\c 2 — 2(3k 3 x 2 + 2k i x + Xi)cd 

+2(k x x + k 2 )ad - 2ab = 0, 

ca t — act — (k 3 x 3 + k A x 2 + Xix)(bc — ad) — (3A; 3 z 2 + 2k^x + \\)c 2 

+(kix + k 2 )ac - a 2 = 0. (48) 

The determining equation for the functions a, b, c and d have now turned out to be nonlinear. 
To solve these equations we further assume that the functions a, b, c and d are polynomials in 
x with coefficients which are arbitrary functions in t. Substituting these forms into Eqs. (|48|) 
we obtain another enlarged set of determining equations for the unknowns and solving the 
latter consistently we obtain nontrivial solutions for the functions a, b, c and d for four sets 
of parametric choices. We present the explicit forms of the associated null function S given 
by ()47|) and the parametric restrictions in Table I. 

Now substituting the forms of S into Eq. and solving the resultant equation we 
obtain the corresponding forms of R. To solve the determining equation for R we again seek 
the same form of ansatz (|28j) but with explicit t dependence on the coefficient functions, that 
is, R = rjn x ^g( t x ^ x y i where Sd is the denominator of S. We report the resultant forms of 
R in Table I. Once S and R are determined then one has to verify the compatibility of this 
set (S, R) with the extra constraint Eq. (j24jl . We find that the forms S and R given in Table 
I do satisfy the extra constraint equation and form a compatible solution. Now substituting 
Si's and i?j's into Eq. one can construct the associated integrals of motion. We report 
the integrals of motion (J) in Table I along with the forms S and R. 

At this stage, we note that the first integral for the case (i) with k 2 ,Xi = has been 
derived in Ref. 18 through Lie symmetry analysis. However, recently, we have derived 9 the 
first integral for arbitrary values of k 2 and Ai. The case (ii) is new to the literature. The 
first integral for the case (Hi) was reported recently in Refs. 9,12 and 13. The first integral 
for the case (iva) is new to the literature. The case r = discussed as (ivb) is nothing but 
the force-free Duffing oscillator whose integrability has been discussed in Refs. 12 and 14. 
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Parametric restrictions, null forms (5), integrating 
x + (kix + fe 2 )x + fe 3 x 3 + k^x 2 + Aix = 



factors (R) and time dependent integrals of motion (I) of 
(identified with the assumed ansatz form of S and R) 



Cases 



Parametric restrictions 



Null form (S) 



Integrating factor (R) 



Integrals of motion (I) 



(i) 



k 2 



kik 2 
3 



(a) I 



(^x 2 -x) 



xe 



33._ 3(-*2T-<) x+klX 2 



3x— — — x+kix 2 



1 3(-fc 2 ±n) 

fc 2 ,Ai t^O, w = - 4Ai)t 



(fei, /c2, Ai : arbitrary) 



(fei, k2, Ai : arbitrary) 



-(fe 2 =F^) + fcia;, 



(fc2±w) . 

e 2 



w = (fe| - 4Ai)2 



m 



fei, k 3 = 0, Ai = ^> 



(fe 2 , &4 : arbitrary) 



(x + &lx) 



iva 



fe 3 = fe4 " ^ 



3 ' 



Ai = 2 -$, r^O 



(fei, fc2, r : arbitrary) 



— + K\X + (3i . +fc2:c) 



(fe 2 x + 3x)e ( - 3 



2(2-r)fc 2 - 



(%a; + rfe 3 x 2 + x) r 



x | x + %x + rfc 3 x 2 



2(2-r) 

e 3 



fe 2 t 



jfe 2 t + log(4fe 2 x + 3feix 2 + 12x) 
4(fe 2 x + 3x) 



(4fe 2 x + 3feix 2 + 12a) ' 



(ivb) 



fei = 0, fe4 = 0, 



(f X + f x + fe 3 x 3 ) 



(x+*f) 



e3 fe2t (x + -|-) 



/ = e3 fc2t 



A 2 , i fe T : _i_ Mt 2 4- k T 4 



Ai = ^,r 







(fe 2 , fe 3 : arbitrary) 



Since we obtained only one integral in each case, (except case (i) where we have found 
second explicit time dependent integral, see Ref. 9), which are also time dependent ones, 
we need to integrate them further to obtain the second integration constant and prove the 
complete integrability of the respective systems, which is indeed a difficult task. 

In this connection we have introduced a new method 1,9 which can be effectively used to 
transform the time dependent integral into a time independent one, for a class of problems, 
so that the latter can be integrated easily. We invoke this procedure here in order to integrate 
the time dependent first integrals and obtain the general solution for all the cases in Table 
I (except case (iv), see below). For the case (iv), we prove the Liouville integrability of it. 

C. Method of transforming time dependent first integral to time independent one 

Let us assume that there exists a first integral for the equation (fTU|) of the form, 



Now let us split the function iq further in terms of two functions such that iq itself is a 
function of the product of the two functions, say, a perfect differentiable function 4(7i(£, x) 
and another function G 2 (t,x,x), that is, 



where iq is a function which involves the variables t, x and x whereas F 2 should involve only 
the variable t and x. We note that while rewriting Eq. (|49|) in the form (|50p. we demand 
that the function F 2 (t, x) in (|4"9*j) automatically to be a function of Gi(t, x). Now identifying 
the function G\ as the new dependent variable and the integral of G 2 over time as the new 
independent variable, that is, 



one indeed obtains an explicit transformation to remove the time dependent part in the first 
integral. We note here that the integration on the right hand side of (joTj) leading to z can 
be performed provided the function G 2 is an exact derivative of t, that is, G 2 = 4;z{t, %) = 
xz x + z t , so that z turns out to be a function t and x alone. In terms of the new variables, 
Eq. ()50|) can be modified to the form 



I = F x (t,x,x) + F 2 (t,x). 



(49) 




(50) 




(51) 




(52) 
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In other words 



Now rewriting Eq. one obtains a separable equation 

£ = /(»>, (54) 

which can lead to the solution after an integration. Now rewriting the solution in terms of 
the original variables one obtains a general solution for the given equation. 

In the following using the above idea we integrate the first integrals given in Table I and 
deduce the second integration constant and general solution. 

D. Application 

Case (ia): k 3 = y, k 4 = ^y 2 -, k\, k 2 and Ai : arbitrary: 

The parametric restrictions given above fix the equation of motion (jTUjl in the form 

x + (hx + k 2 )x + ^-x 3 + + Aix = 0, (55) 

Let us rewrite the first integral associated for this case (vide case (i) in Table I) in the form 



h 



h 2 T^ t 2 



d f , —3 —k 2 ±u. -fc 2 T" 

— ( 1 )e 2 

dt \ k\X 2Ai 



(56) 



where u = \fk 2 — 4A X . Comparing this with the equation (J5U|) . and using (|5*T|). we obtain 

.—3 -k 2 ±u -t 2 T° , , —3 -k 2 ^u. - k 2±«> f 

w=( k7, + ^r )e ' • z=( k7, + ^r )e ' ■ (57) 

Substituting (|57jl into Eq. (|53j). the latter becomes the free particle equation, namely, ^ = 
0, whose general solution is w = Iiz+I 2 , where l\ and I 2 are integration constants. Rewriting 
w and z in terms of x and t one gets 

*(«> = f 6A '' 1 - / f;> | , (58) 

\jfc 1 u>(l + i 1 e wt ) - (A; 2 ± u;) J 2 e ^2-* - fc^l - Ae^) 
where u> = \fk 2 — 4Ai. 
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Interestingly one can consider several sub-cases. In the following we discuss some impor- 
tant ones which are being widely discussed in the current literature. In particular, the dif- 



ference in dynamics arises mainly depending on the sign of the parameter a (= a/A;| — 4Ai). 
We consider the cases (i) k\ < 4Ai (ii) k\ > 4Ai and (Hi) k\ = 4Ai separately. The 
restriction k\ < 4Ai reduces the solution (|58|) to the form 4 , 

x(t) = -, Acos(u t + 5) (5g) 



e 2 



* + 3(4+^g) sin ^ ot + 5 ) ~ ^ cos (<^ + 5))J 



where uiq = — — ^ — ~ an d 5, A are arbitrary constants. A further restriction k 2 = gives us 
the purely sinusoidally oscillating solution 19 



A sin (u; t + 5) . 3a;, 



^)=- — ;rrr u \' ^ o<a<^, w =v%, (eo) 



where A and 5 are arbitrary constants. The associated equation of motion, namely x + 

9 



fcx^x + -^-ar + \\X = 0, admits very interesting nonlinear dynamics, see for example in Ref. 



19. 

On the other hand, in the limit k\ > 4Ai the solution looks like a dissipative/front-like 
one 19 . A further restriction Ai = takes us to the solution of the form 11 

= i jtfg - 1) \ 



Case (ib): k 3 — A; 4 = ■^ 2 -, fcf = 4Ai, fci and fc 2 : arbitrary: 



The third choice k\ = 4Ai in ()58|) leads us to the solution 

, 3(/i+t) y (62) 

Further parametric restriction &2, Ai = provides us the general solution of the form 

*w=r * ( ^'?., y («3) 



The underlying equation, that is, x + fcixi; + ^x 3 = 0, is the I = 1 integrable case of Eq. (JBJ) 
with the solution (JJJ) (see for example in Refs. 18 and 19). 
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Case (ii): k% = 0, k& = ^-(k 2 i — 4Ai), fci, &2 and Ai : arbitrary: 



In this case we have the equation of the form 

fci 



x + (fe a x + k 2 )x + — (k 2 ± \/ fcij - 4Ai)x 2 + Aix = 0. (64) 



The associated first integral reads (vide case (ii) in Table 1) 



i = I x H ^ x + "^"^ 1 e 2 



(65) 



Note that Eq. (|65|) can be rewritten as a Riccati equation of the form 
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. = /e( z^AF^ )t _ ^tVSEJ^j x _ ^. (66) 

The general solution of the Riccati equation is known to be free from movable critical points 
and satisfies the Painleve property. In this sense Eq. (|64|) can be considered as integrable in 
the Painleve criteria sense. However, in the general case, (J66|) . it is not clear whether it can 
be explicitly integrated further. However, for the special case Ai = it can be integrated 



as follows. 



The restriction Ai = ^ fixes the equation of motion (JMj) and the first integral (|65|) in 



9 

the forms 



x + (kix + k 2 )x H x H x = 0, (67) 

3 9 



and 



k 2 k\ .-, \ ■> 

— x H 

3 2 

respectively. Now rewriting in the form f(50j) . we get 



x + ^x+^x 2 e^', (68) 



/ = e^M -^(xe 3 *) + — (xe 3 *)^. (69) 

Identifying the dependent and independent variables from (|6T?|) and using the identities (JSIJ), 
we obtain the transformation 

^2 7- 3 ^2 + 

iw = xe~ , 2; = -— e""3- . (70) 
Using the transformation (J7Uj) the first integral (J68j) can be rewritten in the form 

I = w' + ^w 2 (71) 

21 



which in turn leads to the solution by an integration, that is, 

(21 



w(z) 




tanh 



hi 



(Z - Zq) 



where zq is arbitrary constant. Rewriting (|T2*|) in terms of old variables we get 

3 , / k\I \ , _fc2 



x(t) 
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e ^ 3 )* tanh 



-)(e 



to 



e 3 *) 



where £q is the second integration constant. 



(72) 



(73) 



Case (iii): k±, k 3 = 0, Xi — ^r, k 2 and k A : arbitrary 



The corresponding equation of motion is 



x + k 2 x + k 4 x 2 H -x = 0. 
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(74) 



Rewriting the associated first integral Ji, given in Case (iii) in Table I, in the form (J49|) . we 

get 



1 / 2k 2 \ 6. , kd o 

2 V 5 / 3 



Now splitting the first term in Eq. (J7H|) further in the form (|5()J) . we obtain 



2fc2* f d , 1 2fcat 

I = e 5 —(—=xe 5 
V* V2 



, fc 4 . 

H (xe 

3 V 



(75) 



(76) 



Identifying the dependent and independent variables from (|75jl and using the relations (|51jh 
we obtain the transformation 



1 2fc 2 

u> = —=xe~ 
V2 



5 ^2* 

k 2 



Using this transformation, (|77j) . the first integral (|75|) can be rewritten in the form 



where k\ = 2\f2k i , which inturn leads to 



w ' 2 = Aw A - g 3 , 



(77) 



(78) 



(79) 



where z = 2 A /-j*-z and g 3 



121 

/t'4 



-. The solution of this differential equation can be 



represented in terms of Weierstrass function 12,13 g(z; 0,g 3 ). 
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Case (iv): k 3 = ^ r 2r P k\ , k± = ^y 2 -, Ai = ^jp, fci, &2 and r : arbitrary (but not zero): 



The above parameters fix the equation of motion (fTTH) in the form 

(r- l)k 2 o £4^ n 2H 
x + {hx + k 2 )x + - ' + — - x 2 + — -x = 0, r^O. 

2r z 3 9 

The associated first integral reads (vide case (iva) in Table I) 



(80) 



/ = < 



r ^ 0,2 



(81) 



k 3 



\k 2 t + log(4fc 2 a: + 3^x 2 + 12s) - (4 J!XShU , r = 2 



Rewriting Eq. (JHTj) in the form (|5Uj) . we get 



4r 2 



e 3 



/ = < 



2r 



(xe s**) 2 



r^0,2 (82) 



fe 1 (a;e-r t ) 2 +4^(a : e-r t )e^' 



- log I fc^xe^*) 2 + 4-^(xe^*)e^* 



Identifying the dependent and independent variables from (J82j) and the relations (J51|) . we 
obtain the transformation 

fc 2 j- 3 _ *2 t 



w = xe 3 , 2 = — — e 3 
^2 



(83) 



In terms of the new variables, (JEHJ), the first integral I given above, (J82|l . can be written as 



/ = < 



^k 2 w 4 + w'(w> + ^w 2 ) 



Aw' 



- — log(fciw 2 + Aw') 



, r^0,2 



r = 2. 



(84) 



„ fciui 2 +4u)' 

On the other hand the transformation (}83j) modifies the equation (J80|) to the form 

(r — l)k 2 d 
w" + hww' + v - - J 1 w 3 = 0, r ^ and ' = — . 

2r 2 az 



(85) 
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Finally, for the case r = 0, we have an equation of the form (vide case (ivb) in Table I), 
x + k%x + k%x 3 + \k\x = 0, which is nothing but the force-free Duffing oscillator equation. 
Again using the transformation (|8Hjl . the associated time dependent integral given in Table 
I can be rewritten as 



r W ' 2 i fc3 4 n 
J = — + — w , r = 0. 



(86) 



Though it is difficult to integrate the above time independent first integrals, (|84|). as they 
are in complicated forms, one can easily check that Eq. (|8fi|) (r = 0) can be integrated in 
terms of Jacobian elliptic function 14 and the case r = 1 is already discussed as case (ii) in 
this section. For the other cases one can give a Hamiltonian formulation as in Sec. IIII A 41 
and write the corresponding Hamiltonian as 



where 



H 



(r-1) 



(r-2) 



< >2 p + log( 4^ ; 



k. 2 T 4 1X7 ' 



t / \ l-r 

1 / iU „,2 



, r^0,l,2, 



r = 2 



r = 1 



r = 



p = < 



log(w') 



w 



r = 1 



r = 0. 



(87) 



(88) 



Thus one is ensured of Liouville integrability of system and so (|SU|) for all values of 
r. Further, following the analysis in the above subsection 1111 A 5l one can make a canonical 
transformation (vide Eqs. (}4~2*|) - (pl4*^ to standard nonlinear oscillator equations. 
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E. Summary of results for the q = 1 case: 



To summarize the results obtained in this section, we have identified six integrable cases 
in Eq. (jl(Jj) among which four of them were already known in the literature and the remaining 
two are new. In the following, we tabulate all of them for convenience. 

1. Integrable equations already known in the literature 



(1) x + (kix + k 2 )x + — x 3 + ^^x 2 + X x x = 0, O 

9 3 

(2) x + (k lX + k 2 )x + + ^x = 0, fUD 

3 9 

Qk 2 

(3) x + k 2 x + k 4 x 2 + —^-x = 0, ((71|) 

3 i V J2 



(4) x + k 3 x 6 + k 4 x + Xix = 0. 

We note that the dynamics and certain transformation properties of Eq. (jo3|) have been 
studied in detail by three of the present authors in Refs. 9 and 11 recently. In particular, we 
have shown that this equation admits certain unusual nonlinear dynamics 19 . The dynamics 
of Eqs. (fo7 )) . (f7H) and (133^1 can be found in Ref. 12. 

2. New integrable equations 



(1) x + kixx + k 3 x 3 + Xix = 0, (j33b|) 

bh 2k 2 

(2) x + (hx + k 2 )x + k 3 x 3 + -^x 2 + — = 0, ® 

3 9 

where r 2 k 3 = ^ r and ki,k 2 ,\i and r are arbitrary parameters. We note that (j33b|) 
includes the first equation of MEE hierarchy © as a sub-case. Importantly, we showed that 
(I33b|) is a Hamiltonian system (see Eq. (jHJ)) and so it is Liouville integrable. Equation (|HU|) 
can be transformed to the integrable Eq. ()85j) . Explicit general solution of certain special 
cases, namely, r = 3 or | and r = — 1 or | are reported in Ref. 20. 
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IV. GENERALIZED FORCE FREE DVP FORM OF EQUATIONS 



Let us now consider the case q = 2 in Eq. (jHJ) or equivalently (JUJ), that is, 

x = -{{hx 2 + k 2 )x + k 3 x 5 + k 4 x 3 + X x x) = <f>(x, x). (JTTJ) 

Interestingly Eq. (fTTj) includes another class of physically important nonlinear oscillators. 
For example, choosing k 3 = one can get force-free Duffing-van der Pol oscillator equation. 
With the choice k 2 , k^, Ai = 0, it coincides with the second equation in the MEE hierarchy 
equation. Equation (fTTj) with the restriction fc 3 = k 4 = and Ai = (uJq + ^-), has 
been investigated in a different perspective in Ref. 4. However, the general equation of the 
form has never been considered for integrability test and so we perform the same here. 

To identify integrals of motion and the general solution of Eq. (fTTj) we again seek the PS 
procedure. As the calculations are similar to the q = 1 case of Eq. Q which was carried 
out in the previous section, in the following, we give only the important steps. 



A. The case I t = 

By considering the same arguments given in Sec. IIII A 1\ the null form S can be fixed 
easily in the form 

g_ ((hx 2 + k 2 )x + k 3 x 5 + k 4 x 3 + Aix) , . 

x 

The respective R equation becomes 

R t + xR x — ((k\X 2 + k 2 )x + k 3 x 5 + k 4 x 3 + Xix)R± 

= ((k\x 2 + k 2 ) + ^ l x 2 + k 2 )x + k^ + k i x 3 + X l x) )K m 

To seek a particular form for R one may seek a suitable ansatz. We assume R to be of the 
form (|28|) and investigate the system (|9U|) as before. Following a similar procedure we find 
that a nontrivial particular solution for (|9()j) exists in the form 

R= , 1W , (91) 

(^(^ + k 2 x) + xY 

where r, k\ and k 2 are arbitrary parameters and the remaining parameters, k 3 ,ki and Ai, 
are fixed by the relations 

fc 3 = ^*?, fa=%^fafe, A 1= ^4 (92) 
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Further, we confirmed the compatibility of the functions S and R with the extra constraint 
(124)1 also. We note that unlike the earlier case, q — 1, we do not get a nontrivial solution 
for the parametric restriction k 2 , A; 4 = 0. The above restrictions fix the Eq. (jll)) to the 
following specific forms: 



(ia) x + {kix 2 + k 2 )x + 



(r - 1) ^ + ^r-ljM* ^ + (r - l)k 2 



3r 2 



3r 2 



(ib) x + k%x 5 + k^x 3 + X\x = 0, r = 



2 x = 0, r ^ 0(93a) 
(93b) 



Now substituting ()H9*j) and (j9*T)l into (|23|) and evaluating the integrals we obtain the first 
integrals in the form 



[ia] 



(ib) 



(r — i) ,h o , . 
x + + k 2 x) 



x(x + —x 3 + k 2 x) + ^ 9 ^ (^-x 3 + fc 2 £) 2 



, r^0,2, (94a) 



Ii 



6x 



(6x + 3k 2 x + kix 3 ) 



r 2 3 

log(6i + 3k 2 x + kix 3 ), r = 2, (94b) 



x 



2 6 4 2 ' 



(94c) 



Further, as in the q = 1 case in Sec. IIII A 4) the integrals ()94jl can be recast into the 
Hamiltonian form 



\ia) 
(ib) 
(ic) 
(ii) 



H 



r — l)p 



T-1 

r-1 



(r-2) 



p(—x 3 + k 2 x) 
r 3 



r^0,l,2, 



.6, 



if = — + — x p + log(-), r = 2, 
2 6 p 

ki 

e p + —x 3 + k 2 x, r — 1, 



H = — + — x D + —x 4 + — aT, r = 0. 
2 6 4 2 ' 



where the corresponding canonical momenta respectively are 



(ia, b) 
(ic) 
(ii 



V 



(r-1) 
p = logx, r = 1 



(r - l) f h 
~3 



x + - L (—x 3 + k 2 x) 



(l-r) 



r^0,l, 



(ii) p — x, r = 0. 



(95a) 
(95b) 
(95c) 
(95d) 



(96a) 
(96b) 
(96c) 
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Note that in the above the parameters r, k\, k 2 , k% and Ai are all arbitrary. We also note 
here that unlike the q — 1 case discussed in Sec. IHIl so far we have been unable to find 
suitable canonical transformations for the above Hamiltonian systems so that the standard 
'potential' equation results. The problem is being further investigated. 

B. The case I t ^ 

Now let us study the case I t 7^ 0. In this case S has to be determined from Eq. (J22J), that 



As we did in the q = 1 case of Eq. Q we proceed to solve Eq. (}9T)) with the same form of 
ansatz (jTfjl . Doing so we find that Eq. f)97j) admits non-trivial forms of solutions for certain 
specific parametric restrictions. We report both the parametric values and their respective 
forms of S in Table II. 

Now substituting the forms of S into Eq. (|2^j) and solving the resultant equation we 
obtain the corresponding forms of R. Once S and R are determined then one has to verify 
the compatibility of this solution with the extra constraint ()24|). Then one can substitute 
the null forms and integrating factors into (|23|) and construct the associated integrals of 
motion. We report the integrating factors (R) and time-dependent integrals of motion (J) 



The remaining task is to derive the general solution and establish the complete integra- 
bility of Eq. for each parametric restriction. We again adopt the procedure given in Sec. 
1111 CI and transform the time dependent integrals into time independent ones and integrate 
the latter and deduce the general solution. As the procedure is exactly the same we provide 
only the results in the following. 

Case (ia): k 3 = ^, & 4 = ^p, k\, k 2 and Ai : arbitrary: 

Substituting the parametric restrictions given above in Eq. (fTTj) . we get 



is. 



S t + xS x — {{kix 2 + k 2 )x + k 3 x 5 + k^x 3 + Xix)S± 

= (2k l xx + 5k 3 x 4 + 3k 4 x 2 + A x ) - (k ± x 2 + k 2 )S + S 2 . 



(97) 



in Table II. 



x + (hx 2 + k 2 )x + -^x 5 + 




x A + Aix = 0. 



(98) 
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TABLE II: Parametric restrictions, null forms (S), integrating factors (R) and time dependent integrals of motion (I) of 
x + (kix 2 + fc 2 )x + &i 3 x 5 + /C4X 3 + X\x = (identified with the assumed ansatz form of S and R) 



Cases Parametric restrictions 



Null form (S) 



Integrating factor (R) 



Integrals of motion (/) 



(i) 



&3 — iQ) /C4 — ^ a 



(ki, k2, Ai : arbitrary) 



ki ^3 



x — X 



xe 



( n \ T — „^u>t ( 4x+2(k 2 ±u)x+k 1 x 3 \ 
\ u ) 1 ~ e \4x+2(k 2 T^)x+k 1 x- i ) ' 



(X 



fc 2 ,Ai 7^0, w = (fc|-4Ai)3 



(6)/ 



-t + 



k\ = 4Ai 



(ii) 



fc 3 = 0,fc 4 = ^(fc2±w), 

(fei, fc 2 , Ai : arbitrary) 



■(/c 2 =F w) + fcix 2 



e 2 



k T 3l P (^)t 



w = (k% - 4Ai)2 



^2 ™ I ^2 ™ I , ,-y>3 

X+ ^X 



(iii) 



fei, fc 3 = 0, Ai = ^3- 
(fc 2 , &4 : arbitrary) 



(x + ^x)el fc2i 



I = e-i 



k 2 t 



X± _i_ ^2™™ 1 ,. 



18" 



4 



/ = (ifx 6 + (x + %x)(x + ^x + f X 3 ) j 



(iva) 



(ki, k2, r : arbitrary) 



4 (4x + fc 2 x) 



3(2-r) 



(& 2 x + 4x)e 4 



(^f-x + r/c 3 x 3 + x) r 



c + ^x + r/c 3 x 3 ) e 3(2 4 r)fc2 *, r/2 



X X 



J = |fc 2 t + log(6fc 2 x + 4A;ix 3 + 24x) 



6(/c 2 x + 4x) 



(6fc 2 x + 4fcix 3 + 24x) ' 



r = 2 



(ivb) 



k\ = 0, &4 = 0, 

A 1 = ^,r = 
(& 2 , fc 3 : arbitrary) 



x+ ^x 



e 2 (x + — x) 
4 



We observed that the first integral of this case (i) (see Table II), when rewritten, is nothing 
but the Bernoulli equation which can be integrated strightforwardly 21 and it leads to the 
general solution of the form 

x(t ) = ( sk^-hf y 

KlfhH-h +u)- e^k x k 2 {k 2 + + 8/ 2 A; 2 A 1 e( fe 2+^ + %I x k x \ x &* ) 



where u = \Af — 4Ai- A sub-case of the Eq. ()98|) . namely, k\ < 4Ai has been studied by 
Smith 4,22 who showed that the corresponding equation of motion admits a damped oscillatory 
form of solution, namely, 

x(t) = Acos^t + 5) ^ 



,k 2 t _ kiA , kiA 

4k 2 ^ 4(fc2+4o;2) 



2ujQsin2(uot + 5) — k 2 cos2(u t + 5) 



where ojq = | y4Ai + k\ and 5, A are arbitrary constants. 

On the other hand for k\ > 4Ai, the solution (|9*9*|) becomes dissipative type having a 
front-like structure. In particular, for Ai = we get a solution of the form 

x(t) = ( iVhji^ - 1) \ 

V (-h + 2A; 1 / 1 e fc ^(2 + k 2 I x te k *) + 4k 2 I 2 e 2k ^ J 



Case (ib): k 3 = k^ = ^p, k\ = 4Ai, k\ and k 2 : arbitrary: 



In this case we get the general solution of the form from (jlOlj) as 

x(t) = ( 2 ill+ll ) h (m) 

W \2e k ^I 2 - h(2 + I\k\ + 2k 2 t + k%t 2 + 2I x k 2 {\ + k 2 t))J V ; 

One may note that a sub-case of this equation, namely, k 2 = Ai = leads us to the second 
equation in the MEE hierarchy (JBJ) and the corresponding solution follows from Eq. (jl(J2j) as 

«^(^rj»)' (103) 

This form exactly coincides with the solution (J7J) for / = 2. 



Case (ii): k 3 = 0, k^ = ^-(k 2 ± a/ k\ — 4Ai), ki, k 2 and Ai : arbitrary: 
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The repective equation of motion and the first integral are (see Table II) 



x + (hx 2 + k 2 )x H (A; a ± \ k\ - 4Ai)x 3 + \ x x = 0, (104) 

6 v 



and 



7=laH ^ x + —x* je = 

Eq. (|1U5|) can be rewritten as an Abel's equation in the form 



(105) 



± = /e< s ) f - ( 2 -Jx-yX 3 . (106) 

It is not clear whether Eq. (jlOfijl can be explicitly integrated in general. However, for the 
special case Ai = ^k\ it can be integrated as follows. 

The restriction Ai = ^ fixes the equation of motion ()104|) and the first integral ()1U5|) in 



the forms 



x + (hx 2 + k 2 )x + ^x 3 + ^x = 0, (107) 
4 lo 



and 



x + ^x+^x 6 )e^\ (108) 

respectively. 

Now following our procedure given in Sec. 3.3 one arrives at the general solution 1 as 

1, / (w-a) 2 \ r- f-wy/3 



a 

Z + Zo = -3I 



( (w-a) 2 \ r— 

o lo s + ^ 3 arctan i 

2 \w 2 + aw + ar J \ 2a + w 



(109) 



with w = xe~ , z = —-j^e 2 and a = \ and zo is the second integration constant. 
Rewriting w and z in terms of old variables one can get the explicit solution. 



Case (iii): fci, k$ — 0, Ai = A; 2 and A; 4 : arbitrary: 



The parametric restrictions given above fix the equation of motion (fTTjl to the force-free 
Duffing oscillator, namely, x + k 2 x + k 4 x 3 + ^-x = 0. We have already discussed the general 
solution of this equation in Sec. IIHI (vide case (iv)). 
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Case (iv): k% = - 3 *) fcl , k± = ^p, Ai = ^r, fci, &2 and r : arbitrary: 



The equation of motion turns out to be 



x + (kix + k 2 )x + 



(r — l)k\ 5 fciA^ 3 3fe 



x b + -^ar + —fa; = 0, r ^ 0. 



(110) 



3r 2 ~ ' 4 " 16 

Rewriting the associated first integral /, given in Case (iv) in Table II, in the form (|50p. we 
get 



x | ^(xe 4 c )e 2 + 



3r 



/ = < 



6^(raT') e T f 
fei(o:e^ t ) 3 +6j;(a:e^ t )eTr < 



log(A; 1 (xe" 



xe 4 



6^(xe 4 *)e 2 



r^0,2, 



r = 2 



r = 



(111) 



and identifying the dependent and independent variables from (|111|) and the relations (|51|h 
we obtain the transformation 



^2+ 2 _^2, 

w = xe 4 , 2 = — — e 2 . 



(112) 



In terms of the new variables ()112j) the first integral I given above, (jlllj) can be written as 



w'(w' + ^w 3 ) + %=£klw 6 



L - \og(hw 3 + 6w') 



2 T g ^ , 



, r^0,2 



r = 2, 



r = 0. 



(113) 



On the other hand substituting the transformation ()112|) into the equation of motion 
(HEH we get 



w" + fciwV + — — ^-w 5 = 0, r ^ 0, 
3r 2 



d 

dz 



(114) 
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In the case r = 0, we have an equation of the form (vide case (ivb) in Table II) 



x + k 2 x + k^,x b + 



3^2 

~16~ 



-x 



0. 



(115) 



We note that the Eq. ()114j) is the I = 2 case of Eq. (JEJ). As we mentioned in the introduction 
the general solution of this equation can be obtained only for certain specific choices, namely, 
i. This in turn gives r = 4ki or ~k\. The respective solutions for these values 



3r-' 



of r of Eq. (|114jl can be fixed from Eq. (JJJ) with 1 = 2. The other cases do not seem to 
be amenable to explicit integration. However, all of them can be recast in the Hamiltonian 
form as we see below. 

As the first integrals ()113|) are now 'time' independent ones, one can give a Hamiltonian 
formalism for all the integrals ()113j) by following the ideas given in Sec. IHI A 41 Doing so we 
obtain 



(r-l)p 



(r-2) 



h= | w 3 p + log( p ; 



e p + f w 3 , 



K 2 



r^0,l,2, 



r = 2 



r = 1 



r = 



(116) 



where 



(l-r) 

1 / „,/ | (r-V U .,,3 \ , r ^0,l 



(r-1) \ W> + V fc l W 



logw , 



2 ~ 6 ' 



r = 1 



r = 0. 



(117) 



In this sense these cases may be considered as Liouville integrable systems. Finally, 
for r = case in Eq. f)113j) can be integrated in terms of Jacobian elliptic function (see for 
example in Ref. 23). Again, here, we have not been able to identify canonical transformations 
which can lead to the identification of suitable 'potential' equations. 
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C. Summary of results in q = 2 case: 



To summarize the results obtained for the q = 2 case, we have identified six integrable 
cases in Eq. among which three of them were already known in the literature and the 
remaining three are new. In the following, we tabulate both of them. 

1. Integrable equations already known in the literature 



(1) x + (k lX 2 + k 2 )x + bh. x * + = o, dinzD 

4 16 

2k 2 

(2) x + k 2 x + k 3 x 3 + — ^x = 0, (118) 

9 



(3) x + k 3 x 5 + kiX A + Aix = 0. (J93 

We mention that Eq. (jl(J7p is nothing but the force-free DVP whose integrability is estab- 
lished in Ref. 1 and Eq. ()118|) is nothing but the force-free Duffing oscillator 12 ' 14 . 

2. New integrable equations 



4(r - l)hk 2x3 (r- l)k 2 



x + (k x x 2 + k 2 )x + k 3 x° + — ir x s + x = 0, r^O 



k 2 k k 

(2) x + (hx 2 + k 2 )x + -±x 5 + ^-x 3 + Aix = 0, 

lo 4 

fci k-> 3k 2 

(3) x + (k lX 2 + k 2 )x + k 3 x 5 + -^x 3 + —fx = 0, ffTTOil 

4 16 

where r 2 k 3 = — and hi, k 2) Ai and r are arbitrary parameters. We proved that Eq. (J93a|) 
is Liouville integrable one. As far as equation ((5S}) is concerned we derived the general 
solution for arbitrary values of ki, k 2 and Ai- Finally, for Eq. ()110|) though we identified 
only one time dependent integral, we have demonstrated that it can be transformed into 
time independent Hamiltonian and thereby ensuring its Liouville integrability. 
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V. EXTENDED PRELLE-SINGER METHOD TO GENERALIZED EQ. © 

One can investigate the integrability properties of Eq. Q by considering the cases 
q = 3,4,5, ... , one by one and classify the integrable equations. Since the procedure and 
the mathematical techniques in exploring the integrating factors (R), null forms (S), first 
integrals (/) and general solution are the same in each case we do not consider each case in 
detail. We straightaway move to the case q = arbitrary, that is, q G R and not necessarily 
an integer, and present the outcome of our investigations. 

As we did earlier, we consider the cases I t = and I t ^ separately for the choice 
q = arbitrary also. First let us consider the case I t = 0. 



A. The case I t = 

By considering the same arguments given in Sec. 3.1.1 the null form S and the integrating 
factor R can be fixed easily in the form 

((A^x 9 + k 2 )x + k 3 x 2q+1 + k 4 x 1+q + A x x) 



S 



x 



R = —— — , (119) 

(^(j^x^ + k 2 x) + xY 

where k\ and k 2 are arbitrary and the remaining parameters, ^3,^4 and Ai, are related to 
the parameters k\ and k 2 through the relations 

fc 3 = ^-f^ + l)^i 2 , h = ^^-(q + 2)kk2, Xi = ^^kl (120) 
where k\ = z^n- The above restrictions fix Eq. © to the following specific forms: 

{ia) x+{{q + l)^ + k 2 )x + ^-^[(q + l)k 2 lX 2q+1 

+ {q + 2)k 1 k 2 x q+1 + k\x\ = 0, r ^ (USD 
(ib) x + k 3 x 2q+1 + k 4 x q+1 + \ x x = 0, r = 0. (121) 

Now substituting ()119|) into (|23|) and evaluating the integrals we obtain the first integrals 
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of the form 

(ia) h 



i+ (r — 1) + ^ 



+ kix q+1 + k 2 x) + 



(r-1) 



(hx q+1 + k 2 x) 



iib) h 
(U) h 



x 



k 2 , h +1 , 



(x + f x + fxi +1 ) 

^_ i fc3 x 2 (g+ 1 ) + 

2 2(g+l) (g + 2) 



log(± H x H 



, r = 2, 



, r^0,2, (122a) 
(122b) 



k 



' x q+2 + ^x 2 , r = 0. (122c) 



Further, using the above forms of the first integrals, one can show that the equation of 
motion Q, with the parametric restrictions (jl20}) . can also be derived from the Hamiltonians 



(ta) 

(ib) 
(ic) 
(«) 



H 
H 



r — l)p 



r-2 
r-1 



fr-r 



(r-2) 



p(kix q+1 + k 2 x) 



k 2 , k q+1 , 2{q + i) 



2 1 2 P 
H = e p + hx q+1 + k 2 x, r = l, 

2(9+1) 



H =— + 

2 2(g+l 



-.2" 



+ 



(<7 + l 



-x' 



r^ 0,1,2, (123a) 

(123b) 
(123c) 
(123d) 



where the corresponding canonical momenta respectively are 



(ia, b) p 
(ic) 



1 



x + 



(r-1) 



(r-1) 
p = log re, r = 1 



{k lX q+1 + k 2 x) 



(l-r) 



r^0,l, 



(u) p — x, r = 0. 



(124a) 
(124b) 
(124c) 



With the above Hamiltonian formulation, for the parametric set ()120p . the integrability 
of the associated equation of motion is assured for these parametric cases through Liouville 
theorem. 



B. The case I t ^ 



We use the same ansatz and ideas which we followed for the q = 1 and q = 2 cases to 
determine the forms of S and R. As the procedure is exactly the same as in the earlier cases 
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we present the parametric restrictions and the respective form of expressions of the inte- 
grating factors, null forms and integrals of motions in Table III without further discussion. 



Since we derived only one integral, which is also a time dependent one for each parametric 
restriction, we need to integrate each one of them further and obtain the second integration 
constant in order to prove the complete integrability of each of the cases reported in Table 
III. In the following we deduce the second integral and general solution by utilizing the 
proceduce given in Sec. IHI CI 

Case (ia): k 3 = , , k 4 = ^fy, k 2 and A x : arbitrary: 



We have an equation of the form 

St + {{q + 2)k 1 x q + k 2 )x + k\x 2q+l + k^x^ 1 + \ x x = 0, (fT3J) 

where ki — (q + 2)k x . The corresponding first integral given in Table 3 is nothing but the 
Bernoulli equation which can be solved using the standard method 21 . The general solution 
turns out to be 

*« = (e»'-/ 1 )(eI*-)'(/ 2 + M / g=£)\ft)) \ (125) 

where u = sjk\ — 4Ai. We note here that a sub-case of the above, namely, k\ < 4Ai, has 
been studied by Smith 4 who had shown that the corresponding system admits the general 
solution of the form 

fc 2 , / ~ f -qko . \ « 

x(t) = Acos(to t + 5)e~^' I 1 + qk x A / e— t cos q (uj G t + 5)dt I , (126) 

where uo = | y4Ai + k\ and 6, A are arbitrary constants. For k\ > 4Ai, the solution 
become a dissipative type/front-like structure. In particular, for Ai = the general solution 
takes the form 



x(t) = [ e k2t h 



qk2t > h + hq / h-e- k2t ) dt 



(127) 
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TABLE III: Parametric restrictions, null forms (S), integrating factors (R) and time dependent integrals of motion (7) of 
x + (k\x q + k,2)x + k%x 2q+1 + k^x q+1 + \\x = (identified with the assumed ansatz form of S and R) 



Cases Parametric restrictions 



Null form (S) 



Integrating factor (R) 



Integrals of motion (7) 



2 ^ q+2" 

(-k 2 ±u>) fc! +1 
2 X+ 9 +2 X 



(i) 



(9+2) 



k\k 2 



(a) 7 



= +u;t 



' 9^1 „,o+l _ 
-(9+2) X 



xe 



x 



(x 



(fc 2 ±<+>) 



x + 



ki 

(9+2) 



fc 2 ,Ai t^O, lo = (fc|-4Ai)3 



(fei, ^2, Ai : arbitrary) 



(6)7 



-t + 



(¥^W^) ; 



4Ai 



(") 



/C4 



fcl(fc2±ti>) 
2(9+1) 



, fcs = 0, 



-(fe 2 tw) + hx g , 



e 2 



- i k 2 ^ui j 

X + 2 X + ( ? +l) 

uj = (k\ - 4Ai)2" 



(fei, ^2, Ai : arbitrary) 



,fc 2 ±u 



e \ 2 



2fc 2 A , 4fcfe /., x a+T 



(iii) 



ii n \ 2(g+2)fc| 
fcl,fe 3 = 0, Ai = \ q+ ^y2 2 

(^2) ^4 : arbitrary) 



(x + 



2k 2 x ' 
(9+4) < 



2fe 2 x ^fei 
(x + ^ Iy )e (^) 



2(g+2) 



fc 2 t 



x I a^2^^ _i 

2 "T (q+4) "T (^+iF 



(9+2) 



(iv)a 



, _ (r-l)fc 2 
3 ~~ (q+l)r 2 ' 



fe/1 



k\k 2 
(9+2) : 



fc 2 



(9+2) 



+ feix' i 



2q+l 



(fe 2 x + (g + 2)x)e 



(<j + l)(2-r) 



fc 33 .2( 9 + l) fc2a , ^ fcixg + 1 

(9+1) + i X + ^+2 A X + i+2 + 



(i+^F2)^) 



( (^2)^ + rk 3 x q + 1 + x)' 



, _ (g+i)fc| / n 

A l - ( g +2)2 . r ^ U 

(fei, ^2, f ■ arbitrary) 



U i \ (g + l)(2-r) , . 

^x + rfesx^ 1 +x J e («+ 2 ) fe2 *, r / 2 



/ = g±fe 2 t + log(feix« +1 + 2(g + l)(x + fex)) 
2(g + l)(x + ^_x) 



fcixg+ 1 +2(g+l)(x + ^x) 



k 3 x 2q+i 



;iv)b 



fei =0, fc4 = o, 

Al - (q+2)2 . r - U 

(fc2, k% : arbitrary) 



+ 



(q + 2) (x + 



fc2 ■ 
(9+2)' 



(2 9 +2)fc 2 

e (9+ 2 ) (x + 



fe 2 



(9 + 2) 



x 



7=K + 



k 2 xx 



Tq+2) + 2( 9 "+2)2 -T (2g+2) 



+ 



kzx 2 



2q+2 \ Vi+ 2 ) k 2 t 



e (9+2) 



Case (ib): k% = j|, k^ = ^p, &2 = 4Ai, fci and A;2 : arbitrary: 

A general solution for this case can be fixed from ()127j) as 

x(t) = (Ji + t) e -^* f J a + qk J e~ qJ ^\h + *)«dtj ' . (128) 
On the other hand the general solution for the parametric choice k 2 , \\ = turns out to be 

x(t)= (, (« + !)(/. + «)• (129) 

which exactly coincides with the result (JJJ) obtained by Feix et al. 3 for integer q(— I) values. 
Case (ii): k 3 = 0, k A = 2 (q+i) (^2 ± \A1 — 4Ai), fci, A; 2 and Ai : arbitrary: 
The associated equation of motion and the first integral are (see Table III) 



k\ I 

x + {{q + l)hx q + k 2 )x + -f(k 2 ± Jkj - 4Ai)x 9+1 + X t x = 0, (130) 



and 



i — ix+ h±VSUh x + k lX ^) e (^—) S (131) 



where k\ — (g + l)fci- Like in the earlier cases, that is, q — 1 and q = 2, we are able to 
integrate the first integral (j!31j) explicitly only for a specific parametric restriction, namely, 
■^l = (<? + l)^i> where fc 2 = (g + 2)k 2 . In this case the equation of motion (jl30p and the first 
integral, Eq. (|131jh can be recast in the form 

x + (hx 11 + (q + 2)k 2 )x + k 1 k 2 x q+l + (q + l)kjx = 0, (H2) 

and 

/ = (x + k 2 x + k lX q+l ^j e ^ kt , (132) 
respectively. Now comparing ()132|) with (J3UJ), we get 

/ = e ^ 2 ' (J-(xe k2t ^J + k(xe k2t Y q+1) . (133) 
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Next identifying the dependent and independent variables from (|133|) using the relations 
(|51|). we obtain the transformation 

1 



w = xe k2 \ z 



Using the transformation (jl34j) the first integral (j!33|) can be rewritten in the form 



qk 2 



I = w' + 

which in turn leads to the solution by an integration, that is, 

dw 



z- z 



I - kMi+V 

where Zq is an arbitrary constant. Solving Eq. (jl36J) we get 24 



z - z 



2 



g-1 
2 



q 



2 ^+ 2i 2 ^ 

> K COS 7T H > Qi 

+ 1 ^ + 1 « + 1 

i=0 y y i=0 

1/^(1+™) 



2z 

sin 7r 

g+ 1 



i 



) Q -8, positive odd number, 

2 5; 2z + l 2 X ,2z + l 

> Kicos 7r H > ijsm 

q + o + l o + l ^ 

y i=0 y y i=0 

+-^ln(l + to), q-a positive even number, 



q+1 



TV 



where g = y- and 



1 / 2i 
Pi = -In ( w 2 — 2w cos 7r + 1 I , Qi = arctan 

2 V 9 + 1 

-R, = —ln\w 2 + 2u>cos — — — 7r + 1 I , Ti — arctan 



w — COS -TT7r" 

g+1 

sin -+V7T 

W + COS ^TfTT 



sin ^+t-tt 



2 V 9+1 
Rewriting w and z in terms of old variables one can get the explicit solution. 



(134) 



(135) 



(136) 



(137) 



Case (iii): fci, k% = 0, Ai = 2 |^^2 2 , &2 and &4 : arbitrary 



The parametric choice given above fixes the equation of motion of the form 

x + (q + 4)k 2 x + A; 4 x (9+1) + 2(q + 2)&: 2 x = 0, (H3) 

where A; 2 = (g + A)k 2 . Rewriting the first integral I given in Case (iii) in Table III, in the 
form flU)) , we get 



(138) 
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Now splitting the first term in Eq. 



dt V^2 



further in the form of (J5Uj) 
+ 



(139) 



(9 + 2) \V2 

and identifying the dependent and independent variables from ()139j) using the relations ((EH), 
we obtain the transformation 



x 

w = —^e 
\/2 



2k 2 t 



qk 2 



-qk 2 t 



Using the transformation ()140}) the first integral ()138|) can be brought to the form 

'J±2> 



I = w' 2 + 



2(9 2 )fc 4 , q+2 ) 



(9 + 2; 



(140) 



;i4i) 



Separating the dependent and independent variables and integrating the resultant equation 
we get 

dw 



z- z 



7 - kM q+2) 



(142) 



where k. 



( i+ 2 ) 

1 — 2 ( g+2 ) ^4 an d ^0 is an arbitrary constant. 



Case (iv): k 



{r-l)kj 



[q+l)k 2 2 



1 > K l U . . fclfe2 \ 

3 ( g+ l) r 2 , M (9+2 ) > *1 ( 9+2 )2 ) 



fci, &2 and r : arbitrary: 



The equation of motion in this case turns out to be 

x + {{q + l)kx q + (q + 2)fc 2 )x + (g + 1)( ^ ~ ^ fc 2 z 2g 

r-9 _L P 



+hk 2 x q + k l 2 )x = 0, r^O 



(HHD 



where fci = (g + l)fci, fc 2 = (q + 2)k 2 . Rewriting the associated first integral /, given in Case 
(iv) in Table III, in the form ((50|) . we get 



( f(r-l)kl 



(xe ht f q+l) + ^(xe k2t ) ( ^(xe fe2t )e 9fc2i + ^(xe* 2 *) 



d 



dt 



d 



k 2 t\ qk 2 t 



^k 2 t\q+l \ c qk 2 t 



dt 



at r 



— ( y xe k2t )e qk2t 
log! — \xe 



/ fc 2 tw+l 



+ ^-(xe ht )e qk2t ) 
dt 



2 l 



I| A( xe k 2t\ I p 2 9 fc 2 t _j_ fc 3 ^ p fc 2 t^2(g+l) 



e ^yK 2 i, _|_ 



2(9+1) 



r ^ 0,2 
(143) 

r = 2 
r = 0. 
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Identifying the dependent and independent variables from (j!43|) and the relations (|5ip. we 
obtain the transformation 

1 



w = xe k2 \ z 



qk 2 



(144) 



d 



Subsituting the transformation (|144|) into (|16|) . one obtains 

w" + {q + l)k lW q w' + (q + 1) ^ ~ - fe> 2g+1 = 0, r ^ 0, (145) 

In terms of the new variables (|144j) change the time dependent first integral into time 
independent ones of the form 



w 



1 rr< ^ 



«/(«/ + kw q+l ) + ^A;> 2(9+1) 



/= < 



£ logK + ^^ +1 ), 



w ' 
< 2 



^ 2( g +l) W 



, r ^0,2, 



r = 2, 



r = 0. 



Once again one can deduce the Hamiltonians in the form 



r (r-l)P 



r-2 



(r-1) 



(r-2) 



0,1,2, 



ff= J i^+V + M^ 1 ), r = 2, 



k 2 ^ 2(g+l) W ' ' U ' 



r = 1, 



with 



p = < 



log(w') 



(1-r) 



r = 1 



u/, r = 0, 

and thereby ensuring liouville integrability of Eq. ()16|) . 



(146) 



(147) 



(148) 
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C. Summary of results in q = 



arbitrary case: 



To conclude the integrability of Eq. 0, we have established the fact that the following 



equations, are integrable 

(1) ti + (k 1 x q + {q + 2)k 2 )x + k 1 k 2 x g+1 + {q + l)klx = 0, (H2D 

(2) x + ((q + 2)k 1 x q + k 2 )x + k 2 lX 2q+1 + k^x^ 1 + A x x = 0, (EH) 

(3) x + (q + A)k 2 x + k 4 x q+1 + 2(q + 2)k 2 2 x = 0, (HJ) 

(4) x + ((q + l)hx q + k 2 )x + ^-^-{(q + l)k\x 2q 

+ {q + 2)k 1 k 2 x q + k\)x = 0, r ^ (JUD 

(5) x + ((g+l)A; 1 a; ,? +(g + 2)A; 2 )i: + (g + l)(A;3X 29 + A; 1 A; 2 a; <? + fc 2 2 )x = 0, (USD 



where r 2 k 3 = (r — 1)A; 2 and k\, k 2 , fc 4 , Ai and r are arbitrary parameters (for simplicity we 
have removed hats in k^s, % = 1,2, in Eqs. fll2|) - (fTH|) ) . The significance and newness of the 
equations (fl~2|) - |T6|) are already pointed out in Sec. II Bl 

VI. DISCUSSION AND CONCLUSIONS 

In this paper, we have investigated the integrability properties of Eq. and shown that 
it admits a large class of integrable nonlinear systems. In fact, many classical integrable non- 
linear oscillators can be derived as sub-cases of our results. One of the important outcomes 
of our investigation is that the entire class of Eq. (jHJ) can be derived from a conservative 
Hamiltonian (vide Eq. (j!23|) ) eventhough the system deceptively looks like a dissipative 
equation. 

From our detailed analysis we have shown that Eq. admits both conservative Hamil- 
tonian systems and dissipative systems, depending on the choice of parameters. As far as 
the former is concerned we have deduced the explicit forms of the Hamiltonians for the 
respective equations. In fact, for the case, q — 1, we have constructed suitable canonical 
transformations and transformed the equations into conservative nonlinear oscillator equa- 
tions. However, the canonical transformations for the conservative Hamiltonian systems for 
the cases q — 2, . . . , arbitrary, if at all they exist, still remain to be obtained. Exploring the 
classical dynamics underlying these conservative Hamiltonian systems is also of consider- 
oble interest for further study. As far as dissipative systems are concerned we have not only 
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shown that Eq. Q contains the well known force-free Helmholtz, Duffing and Duffing-van 
der Pol oscillators but also have several integrable generalizations which is another important 
outcome of our investigations. The study of chaotic dynamics of these nonlinear oscillators 
under further perturbations is one of the current topics 22 in the contemporary literature 
in nonlinear dynamics. In principle one can extend such analysis to the above generalized 
equations as well. 

In this paper, we have also not touched the question of linearizability of the integrable 
cases of Eq. (JOJ). In our earlier work, we have shown that the Eq. ()55|) is linearizable to the 
free particle equation, = 0. Of course one can show that this is the only linearizable 
equation in (JHJ) through invertible point transformation 9 ' 11 ' 18 . However, linearizablity of 
other integrable cases through more general transformations still remains to be explored. 

In addition to the above, we have also carried out the Painleve singularity structure 
analysis of Eq. (JUJ) and compared the results obtained through both the methods. The 
details of this will be published elsewhere. 

As we mentioned at the end of Sec. m the crux of the PS procedure lies in finding the 
explicit solutions satisfying all the three determining Eqs. f!22)l - ()24)l . In this paper we have 
considered only certain specific ansatz forms to determine the null forms S, and integrating 
factors R. As a consequence only a specific class of integrable equations have been derived. 
It is not clear, whether these ansatz forms used in this paper exhaust all possible integrable 
cases of Eq. (JHJ). One needs to consider more generalized ansatz forms, and if possible to 
solve Eqs. (j2^|) - (|2^1 for the most general forms of R and S, and try to identify all possible 
integrable cases underlying Eq. (jHJ). This is being explored further. 

VII. ACKNOWLEDGMENTS 

The work of VKC is supported by CSIR in the form of a Senior Research Fellowship. The 
work of MS and ML forms part of a Department of Science and Technology, Government of 
India, sponsored research project. 



44 



1 V. K. Chandrasekar, M. Senthilvelan and M. Lakshmanan, J. Phys. A 37, 
4527 (2004). 

2 S. N. Pandey, M. Senthilvelan and M. Lakshmanan, Classification of Lie point 
symmetries of nonlinear dissipative system of the type x + f(x)x + g(x) = 
preprint, to be submitted for publication. 

3 M. R. Feix, C. Geronimi, L. Cairo, P. G. L. Leach, R. L. Lemmer and S. 
Bouquet, J. Phys. A 30, 7437 (1997). 

4 R. A. Smith, J. London Math. Soc. 36, 33 (1961). 

5 M. Prelle and M. Singer, Trans. Am. Math. Soc. 279, 215 (1983). 

6 Y. K. Man and M. A. H. MacCallum, J. Symbolic Computation 11, 1 (1996). 

7 L. G. S. Duarte, S. E. S. Duarte, A. C. P. da Mota and J. E. F. Skea, J. Phys. 
A 34, 3015 (2001). 

8 L. G. S. Duarte, S. E. S. Duarte and A. C. P. da Mota, J. Phys. A 35, 1001 
(2002); 35, 3899 (2002). 

9 V. K. Chandrasekar, M. Senthilvelan and M. Lakshmanan, Proc. R. Soc. 
London A461, 2451 (2005). 

10 V. K. Chandrasekar, M. Senthilvelan and M. Lakshmanan, J. Nonlinear Math. 
Phys. 12, 184 (2005). 

n V. K. Chandrasekar, M. Senthilvelan and M. Lakshmanan, Chaos, Solitons 
and Fractals 26, 1399 (2005). 

12 M. Lakshmanan and S. Rajasekar, Nonlinear Dynamics: Integrability Chaos 
and Patterns (Springer- Verlag, New York, 2003). 

13 J. A. Almendral and M. A. F. Sanjuan, J. Phys. A 36, 695 (2003). 

14 S. Parthasarathy and M. Lakshmanan, J. Sound and Vib. 137, 523 (1990). 

15 G. W. Bluman and S. C. Anco, Symmetries and Integration Methods for Dif- 
ferential Equations (Springer- Verlag, New York, 2002). 

16 T. C. Bountis, L. B. Drossos, M. Lakshmanan and S. Parthasarathy J. Phys. 
A 26, 6927 (1993). 

17 V. K. Chandrasekar, M. Senthilvelan and M. Lakshmanan, On the complete 
integrability and linearization of nonlinear ordinary differential equations- part 
II: third order equations submitted to Proc. R. Soc. London A. 



45 



18 F. M. Mahomed and P. G. L. Leach, Quaestiones Math. 8, 241 (1985); 12 
121 (1985); P. G. L. Leach, M. R. Feix and S. Bouquet, J. Phys. A 29, 2563 
(1988); P. G. L. Leach, J. Math. Phys. 26, 2510 (1985); R. L. Lemmer and P. 
G. L. Leach, J. Phys. A 26, 5017 (1993). 

19 V. K. Chandrasekar, M. Senthilvelan and M. Lakshmanan, arXiv:nlin. 
SI/0408054 submitted to Phys. Rev. E. 

20 E. L. Ince, Ordinary Differential Equations (Dover, New York, 1956). 

21 G. M. Murphy, Ordinary Differential Equations and Their Solutions (Affiliated 
East-west press, New Delhi, 1969). 

22 D. L. Gonzalez and O. Piro, Phys. Rev. Lett. 50, 870 (1983); Phys. Rev. 
A30, 2788 (1983). 

23 M. Lakshmanan and J. Prabhakaran, 1973 Lettere al Nuovo Cimento 7, 689 
(1973); M. Lakshmanan, Lettere al Nuovo Cimento 8, 743 (1973). 

24 I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series and Products 
(Academic press, London, 1980). 



46 



